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ABSTRACT 

We  present  results  from  a  deep  Chandra  observation  of  Abell  2052.  A2052  is  a  bright,  nearby,  cooling  flow  cluster, 
at  a  redshift  of  z  =  0.035.  Concentric  surface  brightness  discontinuities  are  revealed  in  the  cluster  center,  and  these 
features  are  consistent  with  shocks  driven  by  the  active  galactic  nucleus  (AGN),  both  with  Mach  numbers  of  approxi¬ 
mately  1.2.  The  southern  cavity  in  A2052  now  appears  to  be  split  into  two  cavities  with  the  southernmost  cavity  likely 
representing  a  ghost  bubble  from  earlier  radio  activity.  There  also  appears  to  be  a  ghost  bubble  present  to  the  NW  of 
the  cluster  center.  The  cycle  time  measured  for  the  radio  source  is  t  ~  2  x  107  yr  using  either  the  shock  separation 
or  the  rise  time  of  the  bubbles.  The  energy  deposited  by  the  radio  source,  including  a  combination  of  direct  shock 
heating  and  heating  by  buoyantly  rising  bubbles  inflated  by  the  AGN,  can  offset  the  cooling  in  the  core  of  the  cluster. 

Key  words:  cooling  flows  -  galaxies:  clusters:  general  -  galaxies:  clusters:  individual  (A2052)  -  intergalactic 
medium  -  radio  continuum:  galaxies  -  X-rays:  galaxies:  clusters 


1.  INTRODUCTION 

The  majority  of  cooling  flow  clusters  contain  powerful  radio 
sources  associated  with  central  cD  galaxies.  Initial  evidence  of 
radio  sources  displacing,  and  evacuating  cavities  in,  the  X-ray- 
emitting  intracluster  medium  (ICM)  was  found  with  ROSAT 
observations  of  a  few  sources  including  Perseus  (Bohringer 
et  al.  1993),  Abell  4059  (Huang  &  Sarazin  1998),  and  Abell 
2052  (Rizza  et  al.  2000).  Models  predicted  strong  shocks  sur¬ 
rounding  the  lobes  which  would  appear  spectrally  hard  in 
X-rays  (Heinz  et  al.  1998).  High-resolution  images  from 
Chandra  have  revealed  many  more  cases  of  radio  sources  af¬ 
fecting  the  ICM  by  displacing  it  and  creating  X-ray  deficient 
“cavities,”  or  “bubbles.” 

The  early  Chandra  observations  revealed  that  the  X-ray- 
bright  rims  surrounding  the  radio  sources  in  cooling  flow 
clusters  were  cooler,  rather  than  hotter,  than  the  neighboring 
cluster  gas  (e.g.,  Hydra  A,  Nulsen  et  al.  2002;  Perseus,  Schmidt 
et  al.  2002;  Abell  2052,  Blanton  et  al.  2003).  The  bright  shells 
showed  no  evidence  of  current  strong  shocks.  More  recently, 
shocks  beyond  the  bubble  rims  related  to  radio  source  outbursts 
have  been  found  in  a  few  clusters  including  M87 /Virgo  (Forman 
et  al.  2005),  Hydra  A  (Nulsen  et  al.  2005a),  Hercules  A  (Nulsen 
et  al.  2005b),  and  MS0735.6+7421  (McNamara  et  al.  2005). 
These  shocks  were  fairly  weak  with  Mach  numbers  in  the 
range  of  1.2-1. 7  (see  McNamara  &  Nulsen  2007  for  a  review). 
Ripple  features  in  the  X-ray  surface  brightness  resulting  from 
the  propagation  of  weak  shocks  or  sound  waves  as  seen  in  a 
long  observation  of  Perseus  (Fabian  et  al.  2003,  2006)  may 
also  contribute  to  heating.  Energy  input  from  buoyantly  rising 
bubbles  of  relativistic  plasma  (e.g.,  Churazov  et  al.  2002),  weak 
shocks  (e.g.,  Reynolds  et  al.  2001),  and  the  propagation  of  sound 
or  pressure  waves  are  able  to  offset  cooling. 

Abell  2052  is  a  moderately  rich,  cooling  flow  cluster  at  a 
redshift  of  z  =  0.0348.  A  powerful  radio  source,  3C  317,  is 
hosted  by  the  central  cD  galaxy,  UGC  09799.  Abell  2052  was 


previously  observed  in  the  X-ray  with  Einstein  (White  et  al. 
1997),  ROSAT  (Peres  et  al.  1998,  Rizza  et  al.  2000),  ASCA 
(White  2000),  and  Chandra  (Blanton  et  al.  2001,  2003). 

We  present  a  deep  Chandra  observation  of  Abell  2052, 
combining  the  earlier  Cycle  1  data  with  data  from  Cycle  6.  This 
longer  observation  reveals  probable  shock  features  exterior  to 
the  bubble  rims  that  contribute  to  heating  in  the  cluster  center. 
We  assume  H0  =  70  km  s-1  Mpc-1,  ClM  =  0.3,  and  QA  =  0-2 
(l"  =  0.69  kpc  at  z  =  0.0348)  throughout.  Errors  are  given  at 
the  90%  confidence  level. 

2.  OBSERVATIONS  AND  DATA  REDUCTION 

Abell  2052  was  observed  with  Chandra  using  the  ACIS- 
S  detector  on  2000  September  3  for  a  total  of  36,754  s  and 
on  2006  March  24  for  128,630  s.  The  events  from  the  2000 
data  were  telemetered  in  Faint  mode  and  the  events  from  the 
longer  exposure  were  telemetered  in  VFaint  mode.  The  data 
were  processed  in  the  standard  manner,  using  CIAO  4.0  and 
CALDB  3.4.3.  After  cleaning,  the  total  exposure  remaining  for 
the  two  data  sets  was  162,828  s.  Background  corrections  were 
made  using  the  blank- sky  background  fields. 

3.  X-RAY  IMAGE 

A  slightly  smoothed  (0M9  Gaussian)  image  of  the  central 
region  of  A2052  is  shown  in  Figure  1.  The  overall  emission  is 
fairly  circular  with  a  slight  ellipticity  in  the  north-south  (N-S) 
direction.  Two  separate  surface  brightness  discontinuities  are 
visible  surrounding  the  cluster  center,  and  these  edges  are  seen 
most  clearly  to  the  northeast  (NE)  at  radii  of  approximately  30 
and  45  kpc  from  the  active  galactic  nucleus  (AGN;  apparent  as 
a  bright  point  source  in  the  center  of  the  image).  Depressions  in 
the  surface  brightness  in  the  cluster  center  have  been  previously 
identified  as  bubbles  inflated  by  the  radio  lobes  associated  with 
the  AGN  (Blanton  et  al.  2001,  2003).  To  the  N,  there  is  a  clear 
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Figure  1.  (f/49  Gaussian-smoothed  Chandra  ACIS-S3  image  of  the  central 
3f4x3f4  region  of  Abell  2052.  The  image  has  been  corrected  for  background 
and  exposure.  Several  bubbles  related  to  AGN  outbursts  are  indicated,  as  well 
as  discontinuities  in  the  surface  brightness  surrounding  the  cluster  center.  These 
discontinuities  are  modeled  as  density  jumps  and  are  labeled  “jump  1”  and 
“jump  2.” 
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inner  bubble,  and  with  the  combined  data  set,  an  outer  bubble 
is  seen  to  the  northwest  (NW).  The  S  bubble  now  appears  to 
be  split  into  an  inner  and  outer  bubble,  with  the  outer  bubble  to 
the  southeast  (SE).  These  outer  bubbles  may  be  ghost  bubbles 
arising  from  a  previous  outburst  by  the  AGN.  The  radio  lobe 
emission  from  the  AGN  fills  the  outer  bubbles  and  the  spectral 
index  steepens  in  these  regions  (Zhao  et  al.  1993).  Bright  shells 
surround  the  inner  bubbles,  and  a  filament  extends  from  the 
N  bright  shell  to  the  AGN.  This  was  found  previously  to  be 
associated  with  Ha  emission  (Blanton  et  al.  2001),  and  this 
feature  is  now  revealed  to  extend  much  closer  to  the  AGN. 

An  unsharp-masked  (0.3-10  keV)  image  of  A2052  is  shown 
in  Figure  2  with  a  slightly  larger  field  of  view  than  the 
image  shown  in  Figure  1 .  It  was  created  by  dividing  an  image 
smoothed  with  a  0.98  arcsec  Gaussian  by  an  image  smoothed 
with  a  9.8  arcsec  Gaussian.  The  features  described  above  and 
visible  in  Figure  1  are  also  seen  in  the  unsharp-masked  image, 
sometimes  more  clearly.  The  outer  bubbles  are  revealed  more 
prominently  and  ripple-like  surface  brightness  discontinuities 
are  seen  surrounding  the  cluster  center,  with  two  discontinuities 
visible  most  clearly  to  the  NE  of  the  AGN. 

4.  DENSITY  JUMP  FEATURES 

To  investigate  the  nature  of  the  surface  brightness  disconti¬ 
nuities,  we  have  fitted  a  spherically  symmetric  density  model, 
including  three  power-law  components  with  two  discontinuities 
to  a  NE  sector  of  the  cluster  and  determined  the  jumps  in  den¬ 
sity  that  are  represented  by  these  surface  brightness  changes. 
The  NE  sector  has  an  opening  angle  of  100°  (position  angles 
—2°  to  98°  from  N).  The  fitted  model  is  shown  with  the  emis¬ 
sion  measure  profile  in  Figure  3.  The  emission  measure  was 
calculated  using  projected  spectral  fits,  and  is  fairly  insensitive 
to  temperature.  Assuming  a  constant  temperature,  rather  than 
fitting  for  the  temperature  in  each  region,  changes  the  emission 
measure  values  by  ~5%.  The  best-fitting  model  is  an  excellent 
fit  with  x2  =  21.2  for  22  dof.  A  single  density  jump  model 
does  not  provide  a  good  fit  with  x2  =  116.4  for  25  dof,  nor 
does  a  model  without  a  jump.  With  the  best-fitting  model,  two 
jumps  in  density  are  found,  at  radii  of  45  and  67  arcsec  (31 
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Figure  2.  Unsharp-masked  Chandra  ACIS-S3  image  of  the  central  4.'2x4.'2 
region  of  Abell  2052. 

and  46  kpc)  from  the  AGN,  coincident  with  the  surface  bright¬ 
ness  edges.  The  power-law  indices  for  the  three  components  are 
—0.34,  —1.25,  and  —1.07,  respectively,  going  from  the  center 
of  the  cluster  outward.  The  inner  density  jump  is  a  factor  of 
1.26^0  q53  and  the  outer  density  jump  is  a  factor  of  1.30^  q75. 
These  density  discontinuities  are  consistent  with  shocks  with 
Mach  numbers  of  1.1 8+ 002  f°r  the  31  kpc  feature  and  L20+00°054 
for  the  46  kpc  feature.  The  inner  feature  was  briefly  noted  in 
Blanton  et  al.  (2001),  and  our  Mach  number  is  consistent  with 
the  limit  found  in  the  earlier  paper.  The  very  similar  Mach  num¬ 
bers  for  the  two  shocks  may  indicate  that  the  feedback  cycle 
for  A2052  is  regular,  with  similar  amounts  of  energy  injected 
into  the  cluster  with  each  outburst  of  the  radio  source.  For  each 
shock,  we  would  expect  the  temperature  to  jump  by  a  factor 
of  1.2. 


4.1.  Spectral  Analysis 

In  order  to  examine  the  temperature  structure  near  the 
potential  shock  regions  as  well  as  in  the  cluster  center  in  general, 
we  generated  a  temperature  map.  The  tessellated  temperature 
map  is  shown  in  Figure  4.  Each  bin  was  fitted  using  only  counts 
from  that  area,  so  the  extraction  regions  are  well  defined  and 
do  not  overlap.  The  bins  were  generated  using  the  algorithm 
provided  by  Diehl  &  Statler  (2006),  which  is  a  generalization 
of  Cappellari  &  Copin’s  (2003)  Voronoi  binning  algorithm, 
and  requiring  roughly  1100  net  counts  per  bin  (see  Randall 
et  al.  2008  for  more  details).  Each  bin  was  fitted  in  the 
0. 6-7.0  keV  energy  range  with  an  APEC  thermal  plasma  model 
with  temperature  and  abundance  allowed  to  vary,  and  absorption 
fixed  at  the  Galactic  value  of  Ah  =  2.71  x  1020  cm-2.  The 
color  bar  gives  the  best-fit  temperatures  in  keV.  No  sharp  rise  in 
temperature  is  seen  associated  with  either  of  the  jumps  in  surface 
brightness  corresponding  to  the  shock  regions.  In  general,  the 
cluster  is  coolest  at  the  center  with  the  bright  shells  of  emission 
surrounding  the  bubbles  exhibiting  the  coolest  temperatures, 
kT  «  0.8-1. OkeV. 

We  have  measured  the  projected  temperature  profile  in  the 
same  NE  sector  used  for  the  emission  measure  profile.  The  fits 
were  performed  using  the  APEC  model,  as  described  above.  The 
profile  is  shown  as  the  solid  points  in  the  upper  panel  of  Figure  5. 
As  in  the  temperature  map,  no  large  jumps  in  temperature 
are  seen  inside  the  shocks,  indicated  with  dashed  lines  in 
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Figure  3.  Emission  measure  profile  for  a  NE  sector  of  A2052  with  a  fitted 
spherically  symmetric  density  model  that  includes  three  power-law  components 
(two  discontinuities).  Jumps  in  density  are  found  at  radii  of  45  and  67  arcsec 
(31  and  46  kpc)  from  the  AGN,  as  indicated  by  the  dashed  lines. 
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Figure  4.  Voronoi-Tesselation  temperature  map  of  the  central  region  of  A2052. 
Contours  of  surface  brightness  are  superposed.  A  large  increase  in  temperature 
is  not  seen  behind  the  shock  regions  in  the  NE  in  this  projected  map. 

Figure  5.  Profiles  of  density  and  pressure  are  also  shown  in 
Figure  5.  The  density  was  determined  by  deprojecting  the  sur¬ 
face  brightness  profile,  assuming  that  the  emissivity  is  constant 
in  spherical  shells.  We  use  the  method  described  recently  in 
Wong  et  al.  (2008),  and  also  in  Kriss  et  al.  (1983).  Since  the 
surface  brightness  is  proportional  to  n2Tl/2,  the  density  is  not 
a  very  sensitive  function  of  the  temperature.  We  have  used  our 
projected  spectral  fits  when  calculating  the  conversion  from 
count  rate  to  flux,  which  is  used  in  combination  with  the  surface 
brightness  deprojection  to  determine  the  density.  As  a  com¬ 
parison,  we  have  also  calculated  the  density  directly  from  the 
deprojected  spectral  fits,  using  the  projected  volume  at  each  ra¬ 
dius,  and  find  that  the  density  values  agree  to  within  ~5%  with 
those  calculated  by  deprojecting  the  surface  brightness.  How¬ 
ever,  since  the  errors  in  the  spectral  fits  (both  for  the  tempera¬ 
tures  and  normalizations)  are  much  larger  for  the  deprojected 
spectra,  the  errors  in  the  density  values  are  much  larger  in  this 
case. 

The  temperature  map  in  Figure  4  and  the  temperature  profile 
(solid  points)  in  Figure  5  both  display  projected  temperatures, 
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Figure  5.  Top  panel:  the  projected  (filled  points)  and  deprojected  (open  circles) 
temperature  profiles  for  the  NE  sector  of  the  central  region  of  A2052.  Center 
panel:  density  profile  for  the  NE  sector.  Bottom  panel:  pressure  profile  for  the 
NE  sector  using  the  projected  (filled  points)  and  deprojected  (open  circles) 
temperatures.  The  dashed  lines  in  all  panels  indicate  the  inner  and  outer  jumps 
found  in  the  fit  in  Figure  3. 


and  so  include  the  contribution  of  gas  along  the  line-of- sight 
to  the  gas  at  the  radii  of  interest,  in  the  vicinity  of  the 
shocks.  This  projection  can  make  it  more  difficult  to  detect 
any  temperature  rise  that  may  be  associated  with  the  shocks. 
We  have  therefore  performed  a  spectral  deprojection  for  the  NE 
sector.  The  spectrum  for  the  outermost  annulus  was  fitted  with 
a  model  including  Galactic  absorption  and  an  APEC  model, 
with  temperature  and  abundance  allowed  to  vary.  The  next 
annulus  in  was  then  fitted  with  a  model  including  the  best  fit 
from  the  outermost  annulus  with  the  normalization  scaled  to 
account  for  the  projection  of  this  outer  annulus  onto  the  current 
annulus  being  fitted.  Parameters  from  the  outer  annulus  fit  were 
frozen  and  another  APEC  component  was  added  for  the  current 
annulus.  Spectra  were  fitted  in  this  manner  for  each  annulus, 
including  the  fits  from  all  exterior  annuli  properly  scaled  and 
frozen.  See  Blanton  et  al.  (2003)  for  further  description  of  this 
technique.  The  temperatures  resulting  from  the  deprojection 
are  shown  as  the  open  circles  in  the  upper  panel  of  Figure  5. 
In  addition,  we  have  used  the  deprojected  temperatures  to 
recalculate  the  pressures,  shown  as  the  open  circles  in  the  bottom 
panel  of  Figure  5. 

With  the  spectral  deprojection,  for  the  inner  shock,  the 
best-fitting  preshock  temperature  is  kT  =  2.95+o\78  keV  and 
the  postshock  temperature  is  kT  =  2.97+q332  keV.  For  the 
outer  shock,  these  values  are  kT  =  3. 26+044  keV  and  kT  = 
3.33+044  keV,  respectively.  Therefore,  for  the  inner  shock,  the 
temperatures  are  consistent  with  a  rise  inside  the  shock  by  a 
factor  as  large  as  1.2,  and  for  the  outer  shock  this  value  is  as 
large  as  1.3.  In  both  cases,  then,  the  temperatures  inside  and 
outside  the  shocks  are  consistent  with  what  is  predicted  for 
shocks  with  Mach  number  1.2  (a  temperature  rise  with  a  factor 
of  1.2). 
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5.  DISCUSSION 


While  the  temperatures  we  measure  in  the  density  jump  re¬ 
gions  are  consistent,  within  the  errors,  with  what  is  expected 
with  a  shock  model,  the  best-fitting  temperatures  are  approxi¬ 
mately  constant  across  the  jumps.  If  a  temperature  rise  is  not 
seen  associated  with  a  potential  shock,  the  abrupt  change  in  den¬ 
sity  may  be  due  to  an  isothermal  shock  or  a  cold  front.  Possible 
isothermal  shocks  have  been  seen  in  Perseus  (Fabian  et  al.  2006, 
Mach  1.26)  as  well  as  the  clusters  A2199  and  2 A  0335+096 
(Sanders  &  Fabian  2006,  both  Mach  1.5).  This  may  indicate 
that  thermal  conduction  is  efficient  in  these  regions.  However,  a 
lack  of  an  observed  rise  in  temperature  may  also  result  from  adi¬ 
abatic  expansion  and  orientation  effects  (McNamara  &  Nulsen 
2007).  Cold  fronts  typically  display  a  sharp  change  in  density  in 
one  direction  (see  Markevitch  &  Vikhlinin  2007  for  a  review), 
and  are  not  usually  seen  to  extend  360°  around  the  center  of  a 
cluster.  The  outer  density  jump  feature  is  much  sharper  and  less 
extended  to  the  NE  than  to  the  southwest  and  could  be  associ¬ 
ated  with  a  cold  front  related  to  either  a  merger  or  sloshing  of 
the  cD  galaxy.  The  inner  jump  feature  extends  360°  around  the 
center  of  the  cluster  and  is  probably  more  likely  associated  with 
a  shock.  However,  this  geometry  could  be  possible  if  there  were 
a  cold  front  occurring  along  our  line-of-sight.  The  central  cD 
galaxy  has  a  fairly  large  peculiar  velocity  of  —296  zb  90  km  s-1 
relative  to  the  cluster  mean  (Oegerle  &  Hill  2001),  which  could 
be  related  to  merger  activity  along  our  line-of-sight  possibly 
producing  a  cold  front. 

We  have  calculated  an  average  sound  speed  of  820  km  s-1 
in  the  central  region  of  the  cluster  based  on  a  temperature  of 
kT  =  2.5  keV.  If  we  assume  the  two  density  jump  features 
are  shocks,  we  can  calculate  a  cycle  time  for  the  radio  source 
(where  the  cycle  time  is  the  interval  between  outbursts).  The 
outer  shock  is  at  a  projected  distance  of  46  kpc  from  the  AGN 
and  the  estimated  speed,  based  on  the  Mach  number  above,  is 
980  km  s_1 .  The  inner  shock  is  at  a  projected  distance  of  3 1  kpc 
from  the  AGN,  and  the  estimated  speed  for  this  feature  is  970  km 
s-1 .  If  the  shocks  are  due  to  a  series  of  outbursts  from  the  radio 
source,  their  positions  and  speeds  of  propagation  imply  a  cycle 
time  of  >  1 .5  x  107  yr.  A  similar  cycle  time  was  inferred  from  the 
Perseus  data  from  measuring  the  ripple  separation  (Fabian  et  al. 
2003).  This  time  is  only  slightly  longer  than  the  synchrotron 
radio  age  of  the  lobes  in  Abell  2052  (9  x  106  yr;  Zhao  et  al. 
1993).  We  can  also  estimate  the  cycle  time  using  the  locations  of 
the  ghost  cavities  and  assuming  that  they  rose  buoyantly  at  some 
fraction  of  the  local  sound  speed.  Using  the  SE  ghost  cavity,  we 
estimate  a  cycle  time  of  approximately  4  x  107  yr  (2  x  107  yr) 
if  the  cavities  rose  at  0.5  times  (1  time)  the  local  sound  speed, 
roughly  consistent  with  our  estimate  from  the  shock  features. 

Fitting  a  cooling  flow  model  (including  Galactic  absorption 
and  an  extra  thermal  APEC  model  to  account  for  overlying  gas) 
to  the  cluster  spectra  within  a  radius  of  137"  =  95  kpc  yields  a 
mass  deposition  rate  of  M  =  55  ±4  M©  yr-1.  This  corresponds 
to  a  cooling  luminosity  of  5.4  x  1043  erg  s_1.  Heating  from  the 
central  AGN  can  come  in  the  form  of  shock  heating  as  well 
as  buoyantly  rising  bubbles  inflated  by  the  radio  lobes.  If  the 
density  jump  features  represent  shocks,  we  may  calculate  the 
shock  heating  per  unit  volume  using 
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where  P  is  the  preshock  pressure,  y  =  5/3,  and  In /go  is  the 
time  interval  between  shocks  (McNamara  &  Nulsen  2007).  For 


the  inner  shock,  using  the  preshock  pressure  of  1 .0  x  10-10  dyn 
cm-2,  and  the  expected  pressure  jump  of  a  factor  of  1.55,  we 
find  Us  =  2.8  x  10-27  erg  cm-3  s-1,  or  an  energy  input  rate  of 
1.0  x  1043  erg  s-1  within  the  spherical  volume  interior  to  the 
inner  shock  using  the  cycle  time  calculated  from  the  potential 
shock  feature  separation  of  1 .5  x  107  yr.  Therefore,  energy  input 
from  the  shocks  may  contribute  to  heating,  but  falls  short  by  a 
factor  of  approximately  5  of  offsetting  the  cooling.  We  may 
also  estimate  the  energy  input  from  buoyantly  rising  bubbles. 
Here,  the  heat  input  comes  from  the  work  done  in  inflating  the 
bubbles  and  the  energy  content  of  any  gas  within  the  bubbles, 
E  =  (5/2 )PV  (Churazov  et  al.  2002),  using  y  =  5/3.  Using 
the  pressure  in  the  bubble  rims  of  P  =  1.1  x  10  10  dyn  cm  2, 
computing  the  N  and  S  bubble  volumes  using  a  bubble  radius 
of  10?5  =  7.2  kpc,  and  using  the  cycle  time  inferred  from  the 
outer  bubble  positions  of  4  x  107  yr  (2  x  107  yr),  we  find  an 
energy  input  rate  of  2.0  x  1043  erg  s  1  (4.0  x  1043  erg  s  l). 
If  we  assume  that  the  bubbles  are  filled  with  relativistic  plasma 
(y  =  4/3)  then  E  =  4P  V  (Churazov  et  al.  2002)  and  the  energy 
input  rate  is  slightly  higher,  3.2  x  1043  erg  s  1  (6.4  x  1043  erg 
s-1).  The  combination  of  direct  shock  heating  and  energy  input 
from  buoyantly  rising  bubbles  can  then  offset  the  cooling  of  the 
gas  in  the  center  of  A2052. 
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